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(54) Method and apparatus for analysing an image 

(57) A method of analysing an image to obtain an 
image value. The image comprises a defined array of 
pixel values. The method comprises 



(1) determining the highest pixel value in the image; 

(2) calculating a lower threshold pixel value from 
the highest pixel value determined in step (1) in 



accordance with a predetermined algorithm; and 

(3) obtaining the image value by statistically analys- 
ing the pixel values in the image which lie in a range 
defined by the lower threshold pixel value calcu- 
lated in step (2). 
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Description 

[0001] The present invention relates to a method and 
apparatus for analysing an image. 
[0002] A number of industrial applications require the s 
analysis of an image of a localised feature. In particular, 
it may be desirable to analyze the image to determine 
the average brightness of the localised feature. How- 
ever, the image of the localised feature may be more 
spread out than the original feature itself (for instance to 
due to errors in the imaging process such as optical 
flare, or light from the original feature being reflected 
over adjacent areas). This causes the problem that the 
average brightness is calculated over an area greater 
than the actual area of the original localised feature, and 
as a result the calculated average brightness is errone- 
ously low. 

[0003] This problem has conventionally been 
addressed by reducing the number of pixels which are 
averaged. For instance, only pixels lying above a fixed 
lower threshold pixel value may be averaged. The dis- 
advantage with this approach is that if the image of the 
localised area has low brightness, pixels which should 
be averaged may fall below the lower threshold pixel 
value. This results in a reduction in dynamic range and 
deviation from linearity in the average brightness meas- 
urement. 

[0004] According to a first aspect of the present inven- 
tion, there is provided a method of analysing an image 
to obtain an image value, the image comprising a 
defined array of pixel values, the method comprising 

(1) determining the highest pixel value in the image; 

(2) calculating a lower threshold pixel value from 
the highest pixel value determined in step (1) in 
accordance with a predetermined algorithm; and 

(3) obtaining the image value by statistically analys- 
ing the pixel values in the image which lie in a range 
defined by the lower threshold pixel value calcu- 
lated in step (2). 

[0005] The invention solves the problem outlined 
above by calculating the lower threshold pixel value 
from the highest pixel value. In contrast to the conven- 
tional approach which uses a fixed lower threshold pixel 
value, the invention uses a variable lower threshold pixel 
value which is related to the highest pixel vaJue in the 
image. This results in increased linearity and increased 
dynamic range. 

[0006] Any suitable statistical method may be 
employed in step (3), but preferably the image value is 
obtained in step (3) by averaging the pixel values. The 
average may be a median or mode average but prefera- 
bly the average is a mean average. 
[0007] Typically, the pixel values are related to the 
intensity of radiation from an original image, and in a 
preferred embodiment the method further comprises 
detecting radiation from the original image to generate 



the array of pixel values. The total radiation intensity 
may be detected, or alternatively the intensity in one or 
more selected wavelength ranges may be detected. 
[0008] A further problem can occur if the image con- 
tains erroneously high pixel values (eg. bright spots). 
This can result in the calculated image value being too 
high. This problem can be dealt with by calculating an 
upper pixel threshold value, and/or by removing high 
pixel values before the algorithm is performed, as dis- 
cussed below. 

[0009] Therefore, the method may further comprise 
calculating an upper pixel threshold value from the high- 
est pixel value determined in step (1) in accordance with 
a second predetermined algorithm, wherein step (3) 
comprises statistically analysing the pixel values in the 
image which lie in a range defined by the upper and 
lower pixel threshold values. In some cases, the bright 
spots may be an order of magnitude higher than the 
image values in the rest of the image. Where the image 
is derived from an assay reaction site, this may be 
caused by cosmic rays. This will result in an erroneously 
large highest pixel value being used. Therefore, in this 
case, the method preferably comprises previously defin- 
ing the defined array of pixel values by reviewing an 
original array of pixel values, and removing one or more 
of the highest pixel values in the original array of pixel 
values. This may be done automatically by statistical 
analysis, or by manual inspection. Both methods 
remove the bright spots from the averaging process and 
consequently improve accuracy. 
[001 0] The predetermined algorithm(s) may calculate 
the threshold(s) as any fixed function of the highest pixel 
value. For instance the thresholds may be calculated as 
the square root of the highest pixel value. Alternatively 
the highest pixel value may be input to a look-up table 
which has been previously loaded with a range of 
threshold values. However, preferably the lower and/or 
upper threshold pixel values are calculated as a prede- 
termined percentage of the highest pixel value deter- 
mined in step (1). 

[0011] Typically, the lower threshold pixel value is 
between 50% and 90% of the highest pixel value deter- 
mined in step (1), preferably substantially 80%. 
[0012] Typically, the upper threshold pixel value is 
between 97% and 99% of the highest pixel value deter- 
mined in step (1). preferably substantially 98%. 
[0013] Background noise may be present in the pixel 
values resulting in a reduction in accuracy and dynamic 
range. Preferably, the background is removed from the 
pixel values, for instance by a method of mathematical 
morphology such as erosion, dilation, opening and/or 
thresholding. The background may be removed before 
or after step (1) or step (3). 

[0014] In accordance with a second aspect of the 
present invention, there is provided a method of analys- 
ing an image to obtain a plurality of image values, the 
image comprising an array of pixel values, the method 
comprising dividing the image into a plurality of regions, 
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and obtaining an image value from each region by a 
method according to the first aspect of the invention. 
[0015] The second aspect of the present invention 
enables an image of a plurality of localised features to 
be analyzed. 

[001 6] Typically, the image is an image of a localised 
feature. For instance the localised feature may comprise 
a feature in a satellite image of the earth's surface, or a 
feature in a telescope image. However preferably the 
localised feature comprises a reaction site containing a 
reactive species which reacts with an analyte and which 
has been exposed to a test sample. 
[001 7] In accordance with a third aspect of the present 
invention, there is provided a method of analysing a 
localised reaction site containing a reactive species 
which reacts with an analyte and which has been 
exposed to a test sample, the method comprising 
detecting radiation from the reaction site to generate an 
image comprising an array of pixel values, the radiation 
being indicative of the presence or absence of the ana- 
lyte in the test sample; and analysing the image by a 
method according to the first aspect of the invention, 
thereby obtaining an image value indicative of the pres- 
ence or absence of the analyte in the test sample. 
[0018] The third aspect of the present invention pro- 
vides improved dynamic range and linearity in the anal- 
ysis of an assay reaction site. Typically, the image value 
is an experimental parameter such as a relative light 
unit (RLU) value. 

[0019] In accordance with a fourth aspect of the 
present invention, there is provided a method of analys- 
ing a plurality of localised reaction sites, each contain- 
ing a reactive species which reacts with a respective 
analyte and each having been exposed to a test sam- 
ple, the method comprising detecting radiation from the 
reaction sites to generate an image comprising an array 
of pixel values, the radiation from each reaction site 
being indicative of the presence or absence of a respec- 
tive analyte in the test sample; dividing the image into a 
plurality of image regions each corresponding with a 
respective one of the reaction sites; and analysing each 
image region by a method according to the first aspect 
of the invention, thereby obtaining a plurality of image 
values each being indicative of the presence or absence 
of a respective analyte in the test sample. 
[0020] The fourth aspect of the invention enables a 
multi-analyte assay to be processed quickly and effi- 
ciently. 

[0021 ] Typically, the or each reaction site is located on 
a solid state substrate, and the or each reactive species 
comprises a ligand covalently bonded to the substrate. 
[0022] Typically, the surface of the substrate between 
the reaction sites is inert with respect to analyte. 
[0023] The image may be generated by illuminating 
the or each reaction site and detecting reflected or 
transmitted radiation. However, preferably the radiation 
comprises chemiluminescent, bioluminescent or fluo- 
rescent radiation emitted from the or each reaction site. 



In this case, a filter may be employed so as to only 
detect the narrow wavelength band of interest. 
[0024] In accordance with a fifth aspect of the present 
invention, there is provided a method of performing an 

5 assay, the method comprising exposing one or more 
localised reaction sites to a test sample, each reaction 
site containing a reactive species which reacts with a 
respective analyte; and analysing the or each reaction 
site by a method according to the third or fourth aspects 

10 of the present invention. 

[0025] According to a sixth aspect of the present 
invention, there is provided apparatus for analysing an 
image to obtain an image value, the image comprising a 
defined array of pixel values, the apparatus comprising 

15 

(1) means for determining the highest pixel value in 
the image; 

(2) means for calculating a lower threshold pixel 
value from the highest pixel value in accordance 

20 with a predetermined algorithm; and 

(3) means for obtaining the image value by statisti- 
cally analysing the pixel values in the image which 
lie in a range defined by the lower threshold pixel 
value. 

25 

[0026] According to a seventh aspect of the present 
invention, there is provided apparatus for analysing one 
or more reaction sites each containing a reactive spe- 
cies which reacts with an analyte and which has been 

30 exposed to a test sample, the apparatus comprising 
means for detecting radiation from the or each reaction 
site to generate an image comprising an array of pixel 
values; and means for analysing the image by a method 
according to the first or second aspect of the present 

35 invention, thereby obtaining one or more image values 
each being indicative of the presence or absence of a 
respective analyte in the test sample. 
[0027] According to an eighth aspect of the present 
invention, there is provided apparatus for performing an 

40 assay, the apparatus comprising means for exposing 
one or more localised reaction sites to a test sample, 
each reaction site containing a reactive species which 
reacts with a respective analyte; and apparatus accord- 
ing to the seventh aspect of the present invention for 

45 analysing the or each reaction site. 

[0028] An example of a method and apparatus 
according to the present invention will now be described 
with reference to the accompanying drawings, in which:- 



50 



55 



Figure 1 is a schematic diagram of an integrated 
analyzer system: 

Figure 2 is a cross-section of an integrated analyzer 
system; 

Figures 3-5 are examples of substrates for use in 
the analyzer system; 

Figure 6 is a flow diagram illustrating the main proc- 
ess steps in an assay; 
Figure 7 illustrates erosion; 
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Figure 8 illustrates dilation; 

Figure 9 illustrates opening; 

Figure 1 0 is a graph of calibration curves for each of 

the 12 individual sulphonamides discussed in 

example 1 ; and 

Figures 11-16 are graphs of standard curves for 
FSH LH and PL assays performed on silicon and 
ceramic substrates and image-processed using 
conventional techniques and using a method 
according to the present invention. 



[0029] Figure 1 is a schematic diagram of an inte- 
grated analyzer system for the simultaneous, quantita- 
tive detection of a wide range of analytes in a murti- 
analyte format. The system may be used for clinical/vet- 
erinary diagnosis or drug screening. The system com- 
prises a substrate platform 1 which carries an array of 
substrates 2. The platform 1 is controlled by X-Y trans- 
lation controls 3 which position the platform at each 
stage of the analytical procedure. The translation con- 
trols 3 may be associated with a stepper motor (not 
shown) to achieve a positional accuracy of 10 milome- 
ters. A dispenser unit 4 comprising fluid reservoir 14 
and pump 10 dispenses fluids such as test samples, 
chemiluminescent reagent etc onto the substrates 2. A 
CCD imaging system S comprising a 2-dimens.onal 
CCD camera 6 images the substrates 2 to generate an 
image comprising an array of pixel values which are out- 
put on output line 18. The translation controls 3. dis- 
penser unit 4. and CCD imaging system 5 are controlled 
by a conventional personal computer 7. which runs 
image processing software to process the images 
acquired by the CCD camera 6. 
[0030] Figure 2 illustrates a system of the type shown 

in the schematic diagram of Figure 1. CCD camera 6 
has a lens 8 which acquires an image of a substrate 2 in 
an imaging position 9. Liquids are dispensed with dis- 
penser pump 10 onto a substrate 2 in a dispensing posi- 
tion 1 1 . The lens 8 and imaging position 9 are enclosed 
inadarkbox12. 

[0031] Substrates 2 of the form illustrated in Figures 
3-5 are typically employed in the system of Figures 1 
and 2 The substrates shown in Figures 3-5 are silicon 
or ceramic chips which have a plane surface which has 
been modified to incorporate a series of d.m- 
ples/wells/p'rts. channels or chambers on the surface. 
These help to increase the surface area of the reaction 
sites, thereby providing the possibility of accelerating 
assay incubation time 

[0032] Referring to Figure 3. a chip 20 has a central 
test sample reservoir 21 to which a liquid test sample is 
added when the chip 20 is at the dispensing position 1 1 . 
Twelve channels 22 radiate from the reservoir 21 to 
respective reaction sites 23. 

[0033] Referring to Figure 4. chip 25 has four test 
sample reservoirs 26-29 to which liquid test sample is 
added at the dispensing position 1 1 . a delivery manifold 
30 test sample delivery channels 31 which deliver the 



test sample to reaction sites (not shown), and a waste 
reservoir 32. 

[0034] Referring to Figure 5. a chip 40 has twelve cir- 
cular reaction sites 41. The chip 40 is a 1cm 2 ceram.c 
5 substrate. .. . 

[0035] Each reaction site 23,41 etc. bears a l.gand 
covalently bonded to the substrate, and each ligand 
reacts with a respective different analyte. The Hgands 
are immobilised in the spatially defined reaction sites by 
w means of microfluidic dispensing of the ligand onto the 
substrate, which is chemically activated. The surface of 
the substrate between the reaction sites is inert with 
respect to the analyte. Suitable methods of preparing 
the reaction sites are described in more detail in co- 
ts pending European Patent Application No. 97302707 >A . 
Covalent immobilisation of the biological ligand in the 
reaction site 23.41 etc is preferred since passive 
absorption interactions are susceptible to changes in 
pH temperature and ionic strengths, and may in some 
20 instances result in release of weakly-bound lisandsdur- 
ing incubation and washing steps, thus contributing to 
poor assay reproducibility. 

[0036] A method of performing an assay using the 
system of Figures 1 and 2 will now be described w* 
25 reference to Figure 6. In a first step 60. a substrate 2 is 
placed on the sample tray 1 . The substrate 2 has a plu- 
rality of reaction sites, each bearing a l.gand which 
reacts with a respective analyte A.B.C.D etc. At step 61 
the sample tray 60 is moved by translation controls 3 
so until the substrate is located in dispensing position 11 
At 62 the dispenser pump 10 dispenses a liquid test 
sample onto the substrate - for instance in the case of 
the substrate of Figure 3 the test sample is dispensed 
into the reservoir 21 . At 62'. the dispenser pump 1 0 dis- 
35 pensesaliouidreagerrtortomesubs^ 

reagent comprises a mixture of labelled reagents 
A B 1 C D" etc which each react with a respective one ot 
the analytes, i.e. A' reacts with A, B" reacts with B eto 
The labelled reagents A'.B'.C'.D' etc also all react with 
40 chemiluminescence reagent. 

[0037] The test sample and labelled reagents then 
pass along the channels 21 into the reaction sites 23, 
exposing the various ligands to the test sample. The 
flow of liquid test sample and labelled reagents may be 
45 achieved by capillary attraction, centrifugal force vac- 
uum force or electroosmoticflow. The use of electroos- 
motic flow would avoid the need for valves, so that no 
moving mechanical parts are used. At 63 the substrate 
is washed to remove excess test sample which has not 
so reacted at a reaction site 23. 

[0038] At step 64 the dispenser pump 1 0 dispenses 
chemiluminescence reagent onto the substrate. This 
reacts with the labelled reagents A.B.C.D etc in a 
chemiluminescent reaction. At 65. the translation con- 
55 trols 3 move the sample tray 1 until the substrates are in 
imaging position 9. At this stage, where an analyte has 
reacted with a ligand at one of the reaction srtes. chemi- 
luminescent radiation is being emitted from that reaction 
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site by a respective labelled reagent. At 66 t the CCD 
camera 6 with a suitable filter is exposed to generate an 
image of the chemiluminescent radiation emitted by the 
chip 2 in the imaging position 9, the image comprising 
an array of pixel values. The array of pixel values is input 
to the personal computer 7 via image output line 1 B, and 
stored in a memory (not shown) for analysis. At 67, the 
personal computer 7, under control of image processing 
software, analyses the image from the CCD camera 6 to 
generate a relative light unit (RLU) value for each reac- 
tion site. At 68, the personal computer 7 outputs the 
RLU values on a monitor or printer (not shown) and at 
69 the personal computer 7 calculates a corresponding 
concentration for each analyte A.B.C.D etc in the test 
sample and outputs the calculated concentration on the 
printer or monitor. 

[0039] The image processing algorithm 67 proceeds 
as follows: 

a) The image is divided into a plurality of regions, 
each region containing one of the reaction sites and 
being defined by a two-dimensional array of pixels. 
The image is divided with reference to a previously 
stored map corresponding with the pattern of the 
reaction sites. 

b) Within each region the computer 7 scans through 
and locates the most intense pixel. Let this intensity 
= X. 

c) Next, it calculates 80% of this pixel intensity 
(=0.8X) and, hence, the top 20% range. 

d) Finally, it identifies all pixels within that area of 
interest that fall within this range and averages the 
intensity of these pixels which it outputs as the top 
20% signal intensity. Pixels that fall outside the 
range are discarded. 

n 

TOP 20% RLU signal ■ 

for all pixel intensities Pj, where 0.8X£Pj*X 
X=top signal, n=no. of pixels within the 80-100% 



range. 

[0040] Other ranges can be used. 40% gives reason- 
ably high intensities. Also, to remove any uncharacteris- 
tically high pixels, ranges of 80%r->98% can be 
examined. 

[0041 ] In an alternative method, in order to deal with 
abnormally high highest pixel values (which may result 
from bright spots caused by cosmic radiation), an addi- 
tional step may be carried out between step a) and step 
b). In this additional step, each region is reviewed to 
determine whether there are any "rogue" values. For 
instance, with an array of one hundred pixel values, 
ninety-nine may lie between 0-100 RLUs and one 
"rogue" pixel may have a value of 2000 RLUs. The 
rogue pixel (or pixels) is discarded before step b). This 



step may be carried out manually, or by a suitable statis- 
tical algorithm. 

[0042] If there is any inherent background in the 
image, this will be included in this signal. If all images 
5 are processed by this method, this will not necessarily 
be a problem but could be restrictive on the dynamic 
range. Therefore, mathematical morphology processing 
is used to remove background. Mathematical morphol- 
ogy is based on the use of set operators (intersection, 
10 union, inclusion, complement) to transform an image. 
The processed image has usually fewer details, imply- 
ing a loss of information, but its main characteristics are 
still present. Once the image is simplified, quantitative 
analysis can be performed. 
75 [0043] The principle of a basic operation transforma- 
tion is based on the choice of a "structuring element", 
characterized by a shape, a size and the location of its 
centre. Each object of the image will be compared with 
this element and this is achieved by moving the element 
20 such that its centre hits all the points of the image. For 
each position of the element, one looks for the inclusion 
of the element with the objects of the image. 

removes isolated points and small parti- 
cles, shrinks other particles, discards 
peaks on the boundary of the objects and 
disconnects some particles. In practice, 
with the structuring element, each central 
pixel is set to the lowest value in that par- 
30 ticular position. Referring to Figure 7, a 

pixel value in pixel location 70 at the cen- 
. tre of a 3x3 structuring element 71 is set 
to the lowest value within the structuring 
element 71 , in this case value "1 
35 Eilafion - fiUs small holes inside particles and gulfs 
on the boundary of objects, enlarges the 
size of the particles and may connect 
neighbouring particles. In practice, within 
the structuring element, each central 
40 pixel is set to the highest value in that 

particular position. This is illustrated in 
Figure 8. 

Opening - is a transformation which does not 
destroy as much information as erosion, 
45 and is the combination of an erosion fol- 

lowed by a dilation. 



[0044] Opening gives an accurate determination of 
the background and can then be subtracted from the 
so original image. Further use of thresholding then pro- 
duces "clean" spots. One iteration is illustrate in Figure 
9. Increased iterations maintain high spot signal and 
give smoother edges. 



55 Thresholding • 



sets the minimum intensity from 
which signal is read. This also helps 
to clean up the image and eliminate 
misleading background. Each spot is 
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individually labelled and the intensity 
from the whole spot averaged. 

[0045] The advantages of morphology are the 
removal of sloping background and the ability to distin- 5 
guish between real signal and surrounding background. 
However, since the whole spot intensity is averaged, the 
actual signal tends to be very low. 
[0046] A way around this problem is to use the mor- 
phology method, but instead of averaging the whole 10 
spot average only the top 20% as outlined previously. 
This method results in high signals, but also an 
improved dynamic range due to the removal of back- 
ground. 

[0047] e.g. consider a top standard of 20.000RLU and is 
a bottom standard of 1200RLU on a background of 
600RLU. 

[0048] With background, dynamic range = 



20000 
1200 



= 167 



[0049] Without background, dynamic range = 



20000-600 
1200-600 ' 



32.3 



[0050] Two illustrative methods of preparing sub- 
strates 2 and performing assays using the system of 
Figures 1 and 2 will now be described below. 

Example 1 Sulphonamide Multi-Analyte Assay 



20 



25 



30 



35 



[0051] In this example 12 individual antibodies, each 
antibody specific for a single sulphonamide, were immo- 
bilised by covalent attachment by contact interactions 
on to discrete regions of a ceramic (aluminium oxide) 
substrate chemically modified epoxide surface. *o 
[0052] The ceramic substances (1cm x 1cm) were 
ultrasonically cleaned using an alkaline detergent 
(RBS35 @ 5% V/V) followed by double deionised water 
and then placed in 6M HCl for 16 hours. The chips 
(1 cm 2 ) were then placed in chromic acid for 1 hour in an 45 
ultrasonic bath. 

[0053] The chips were washed exhaustively with dou- 
ble deionised water and acetone and then dried in an 
oven at 120°C for 2 hours. Following this pretreatment, 
the chips were silanated using the organosilane y-glyci- so 
doxypropyl trimethoxy silane (10%V/V) in anhydrous tol- 
uene, 4-dimethyl aminopyridine (1.25g/L) and 
triethylamine (1% V/V). This mixture was refluxed for 4 
hours and then left overnight at room temperature. The 
chips were washed with toluene and acetone before ss 
curing for 4 hours at 120°C. 

[0054] Following the curing step the chips were placed 
in containers and stored at room temperature until' 



required for spotting of sulphomanide antibodies. The 
sulphonamide antibodies were spotted using a BIODOT 
XY3000 dispenser in the pattern illustrated in Figure 5. 
The sulphonamides assayed were sulphadoxine, sul- 
phamethizole, sulphachlorophyridazine, sulphamethox- 
ypyridazine, sulphamerazine, sulphapyridine, 
sulphisoxazole, sulphathiazole, sulphamethazine, sul- 
phaquinoxaline, sulphadimethoxine, and sulphadiazine. 
[0055] Dispense volumes of approx. 20nl were 
employed for each sulphonamide antibody. The 12 sul- 
phonamide antibodies which formed 12 discrete areas 
on the 1cm 2 chip were incubated for 2 hours at 37°C. 
The chips were washed with phosphate buffered saline 
(pH7.2) containing 2% casein (W/V) and then blocked in 
same buffer overnight at +2-8°C. After washing with 
phosphate buffered saline (PBS) containing PEG300 
(0.05% VA/) the chips were placed in the sample holder 
1. 

[0056] Multi-sulphonamide standards (200^1) and a 
cocktail of sulphonamide horseradish peroxidase conju- 
gates (100ul) were added to the wells of the device as 
appropriate and incubated for 15 minutes at room tem- 
perature. The standards contained 5ng/ml, 10ng/ml, 
50ng/ml and 100ng/ml for each of the 12 sulphona- 
mides. 

[0057] Thereafter the multi-sulphonamide chips were 
washed with PBS/PEG buffer to remove excess rea- 
gents and 300nl chemiluminescent substrate [lumi- 
nol(1.4mM)/urea hydrogen peroxide (9.6mM)] per chip 
was introduced. The chips were imaged using a CCD 
camera device with an exposure time of up to 4 minutes. 
A back-illuminated CCD camera is preferred, to improve 
the capture eff iciency at the wavelength of the light gen- 
erated by the chemiluminescent light reaction (approxi- 
mately 433-445 nm in the case of luminol). The 
standard curves were obtained for each of the 1 2 indi- 
vidual sulphonamides using the new image-processing 
method. The calibration curves for each of the 12 indi- 
vidual sulphonamides are represented graphically in 
Figure 13. The % B/Bo value plotted on the Y-axis rep- 
resents the % inhibition of the zero standard RLU(Rela- 
tive light unit) value caused by each individual 
sulphonamide standard (plotted on the X-axis as log 10 ). 
[0058] This example represents a mutti-analyte assay 
for analytes of low molecular weight, where the assay is 
based on a competitive immunoassay format. 

Example 2 Multi-Analyte Assay For Prolactin (PL), Folli- 
cle Stimulating Hormone (FSH) And Luteinizing Hor- 
mone (LH). 

[0059] In this example a multi-analyte assay was per- 
formed for 3 hormones of large molecular weight This 
example represents a multi-analyte assay for a sand- 
wich based immunoassay. No significant cross-reactiv- 
ity were observed when the three hormones were 
determined in the same panel. The chemical pretreat- 
ment and silanation procedures were exactly as 
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described in example 1. 

[0060] The individual PL.FSH or LH monoclonal anti- 
bodies (approx.20nl antibody dispensed) were immobi- 
lised on discrete areas of the chemically modified 
substrate. The multi-analyte assays were performed on 5 
both silicon and ceramic substrates with an epoxide sur- 
face as described in example 1. 
[0061] In the assay 150uJ of a multiple LH/PL7FSH 
serum based standard and 150uJ of a diluent assay 
buffer were added to the chip and incubated for 1 5 min- 10 
utes at room temperature. Following a wash step 300^1 
of a single conjugate cocktail of LH-HRPO/PL- 
HRPO/FSH-HRPO conjugates was added and incu- 
bated for 1 5 minutes. Thereafter the chips were washed 
to remove excess reagents and the chemiluminescent is 
reagent [luminol(1 .4mM)/urea hydrogen perox- 
ide(9.6mM)] was introduced. 

[0062] The chips were imaged using a CCD camera 
device with an exposure time of up to 4 minutes. The 
standard curves for each of the hormones were plotted 20 
after the images were processed using both the tradi- 
tional morphology/threshold method and the new 
image-processing method. 

[0063] The standard curves are illustrated graphically 
in Figures 10-15. The standard curves obtained using 25 
the traditional morphology/threshold method are 
labelled 80-85. The standard curves obtained using an 
example of an image-processing method according to 
the present invention are labelled 86-91 . It can be seen 
that the linearity and dynamic range indicated by the 30 
curves 86-91 is greiater than the linearity and dynamic 
range indicated by the curves 80-85. The greater linear- 
ity and dynamic range offers a significant time saving to 
the end-user, since a wide dynamic range diminishes 
the need for the retesting of very high concentration test 35 
samples after additional dilution steps, to obtain an 
accurate result. The improved linearity/dynamic range 
also improves intra and inter-assay precision for test 
samples, since a small RLU change does not adversely 
affect results. The increased dynamic range has added 40 
benefits in terms of assay sensitivity, with improved dis- 
crimination of reaction site RLU values, enabling lower 
levels of analyte to be detected. 



Claims 

1 . A method of analysing an image to obtain an image 
value, the image comprising a defined array of pixel 
values, the method comprising 
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(1) determining the highest pixel value in the 
image; 

(2) calculating a lower threshold pixel value 
from the highest pixel value determined in step 
(1) in. accordance with a predetermined aigo- ss 
rithm; and 

(3) obtaining the image value by statistically 
analysing the pixel values in the image which 



lie in a range defined by the lower threshold 
pixel value calculated in step (2) . 

2. A method according to claim 1 , further comprising 
calculating an upper pixel threshold value from the 
highest pixel value determined in step (1) in accord- 
ance with a second predetermined algorithm, 
wherein step (3) comprises statistically analysing 
the pixel values in the image which lie in a range 
defined by the upper and lower pixel threshold val- 
ues. 

3. A method according to claim 1 or 2, wherein the 
lower and/or upper threshold pixel values are calcu- 
lated as a predetermined percentage of the highest 
pixel value determined in step (1). 

4. A method according to claim 3, wherein the lower 
threshold pixel value is between 50% and 90% of 
the highest pixel value determined in step (1). 

5. A method according to claim 4, wherein the lower 
threshold pixel value is substantially 80% of the 
highest pixel value determined in step (1). 

6. A method according to any of claims 3 to 5, wherein 
the upper threshold pixel value is between 97% and 
99% of the highest pixel value determined in step 
0). 

7. A method according to claim 6, wherein the upper 
threshold pixel value is substantially 98% of the 
highest pixel value determined in step (1). 

8. A method according to any of the preceding claims, 
further comprising removing background from the 
pixel values. 

9. A method according to claim 8, wherein the back- 
ground is removed by a process of erosion, dilation, 
opening and/or thresholding. 

10. A method of analysing an image to obtain a plurality 
of image values, the image comprising an array of 
pixel values, the method comprising dividing the 
image into a plurality of regions, and obtaining an 
image value from each region by a method accord- 
ing to any of the preceding claims. 

1 1 . A method according to any of the preceding claims, 
further comprising detecting radiation from an orig- 
inal image to generate the array of pixel values. 

1 2. A method according to any of the preceding claims, 
further comprising previously defining the defined 
array of pixel values by reviewing an original array 
of pixel values, and removing one or more of the 
highest pixel values in the original array of pixel val- 
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ues. 

1 3. A method according to any of the preceding claims, 
wherein the image value is obtained in step (3) by 
calculating the average of the pixel values. 

14. A method according to any of the preceding claims, 
wherein the image has been obtained from a reac- 
tion site containing a reactive species which reacts 
with an analyte and which has been exposed to a 
test sample. 

15. A method of analysing a localised reaction site con- 
taining a reactive species which reacts with an ana- 
lyte and which has been exposed to a test sample, 
the method comprising detecting radiation from the 
reaction site to generate an image comprising an 
array of pixel values, the radiation being indicative 
of the presence or absence of the analyte in the test 
sample; and analysing the image by a method 
according to any of the preceding claims, thereby 
obtaining an image value indicative of the presence 
or absence of the analyte in the test sample. 

16. A method of analysing a plurality of localised reac- 
tion sites, each containing a reactive species which 
reacts with a respective analyte and each having 
been exposed to a test sample, the method com- 
prising detecting radiation from the reaction sites to 
generate an image comprising an array of pixel val- 
ues, the radiation from each, reaction site being 
indicative of the presence or absence of a respec- 
tive analyte in the test sample; dividing the image 
into a plurality of image regions each corresponding 
with a respective one of the reaction sites; and ana- 
lysing each image region by a method according to 
any of claims 1 to 1 4, thereby obtaining a plurality of 
image values each being indicative of the presence 
or absence of a respective analyte in the test sam- 
ple. 

17. A method according to claim 15 or 16, wherein the 
or each reaction site is located on a solid state sub- 
strate, and wherein the or each reactive species 
comprises a iigand covalerrtly bonded to the sub- 
strate. 

18. A method according to claim 16 and 17, wherein 
the surface of the substrate between the reaction 
sites is inert with respect to analyte. 

19. A method according to any of claims 15 to 18, 
wherein the radiation comprises chemiluminescent 
bioluminescent or fluorescent radiation emitted 
from the or each reaction site. 

20. A method of performing an assay, the method com- 
prising exposing one or more localised reaction 



sites to a test sample, each reaction site containing 
a reactive species which reacts with a respective 
analyte; and analysing the or each reaction site by 
a method according to any of claims 1 5 to 1 9. 

> 

21. Apparatus for analysing an image to obtain an 
image value, the image comprising a defined array 
of pixel values, the apparatus comprising 

w (1) means for determining the highest pixel 

value in the image; 

(2) means for calculating a lower threshold 
pixel value from the highest pixel value in 
accordance with a predetermined algorithm; 

is and 

(3) means for obtaining the image value by sta- 
tistically analysing the pixel values in the image 
which lie in a range defined by the lower thresh- 
old pixel value. 

20 

22. Apparatus for analysing one or more reaction sites 
each containing a reactive species which reacts 
with an analyte and which has been exposed to a 
test sample, the apparatus comprising means for 

25 detecting radiation from the or each reaction site to 
generate an image comprising an array of pixel val- 
ues, the radiation from the or each reaction site 
being indicative of the presence or absence of a 
respective analyte in the test sample; and means 
30 for analysing the image by a method according to 
any of claims 1 to 1 4, thereby obtaining one or more 
image values each being indicative of the presence 
or absence of a respective analyte in the test sarrv 
pie. 

35 

23. Apparatus for performing an assay, the apparatus 
comprising means for exposing one or more local- 
ised reaction sites to a test sample, each reaction 
site containing a reactive species which reacts with 

40 a respective analyte; and apparatus according to 
claim 23 for analysing the or each reaction site. 
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Fig.2. 
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Fig.6. 
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Fig.9. 
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Fig.11. 
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Fig.12. 
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Fig. 13. 
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Fig. 14. 
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Fig.15. 

Multianalyte assay Silicon GOPS - Prolactin 

o Morphology Processing (mt=300) 
A New image processing method 



40- 



w 

"D 
C 
CO 
C/3 
3 
O 

sz 

<D 
3 

CO 
> 

—J 



20- 



10 




100 



T 



1000 
Log 10 Prolactin 



10000 



20 



EP 0 902 394 A1 



80- 
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